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Abstract 

We study travelling kinks in the spatial discretizations of the nonlinear Klein-Gordon equa- 
tion, which include the discrete cf> 4 lattice and the discrete sine-Gordon lattice. The differential 
advance-delay equation for travelling kinks is reduced to the normal form, a scalar fourth-order 
differential equation, near the quadruple zero eigenvalue. We show numerically non-existence 
of monotonic kinks (heteroclinic orbits between adjacent equilibrium points) in the fourth-order 
equation. Making generic assumptions on the reduced fourth-order equation, we prove the per- 
sistence of bounded solutions (heteroclinic connections between periodic solutions near adjacent 
equilibrium points) in the full differential advanced-delay equation with the technique of center 
manifold reduction. Existence and persistence of multiple kinks in the discrete sine-Gordon equa- 
tion are discussed in connection to recent numerical results of | ACR03 1 and results of our normal 
form analysis. 



1 Introduction 

Spatial discretizations of the nonlinear partial differential equations represent discrete dynamical sys- 
tems, which are equivalent to chains of coupled nonlinear oscillators or discrete nonlinear lattices. 
Motivated by various physical applications and recent advances in mathematical analysis of discrete 
lattices, we consider the discrete Klein-Gordon equation in the form: 

u n+ i - 2u n + u n -i 

U n = Jo, \- f{U n -l,U n ,U n+1 ), (1.1) 

where u n {t) £ R, n £ Z, t S R, h is the lattice step size, and f(u n -\,u n ,u n+ i) is the nonlinearity 
function. The discrete lattice ( ll.lt is a discretization of the continuous Klein-Gordon equation, which 
emerges in the singular limit h — ► 0: 

utt = Uxx + F(u), (1.2) 
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where u(x,t) G R, x G R, and i G R. In particular, we study two versions of the Klein-Gordon 
equation dl.2t . namely the </> 4 model 

u t t = u xx + it(l - -u 2 ) (1.3) 

and the sine-Gordon equation 

u u = u xx + sm(u). (1.4) 

We assume that the spatial discretization f(u n -x,u n ,u n+ i) of the nonlinearity function F(u) is sym- 
metric, 

f(u n -i,u n ,u n+ i) = f(u n+ i,u n , u n -i), (1.5) 
and consistent with the continuous limit, 

f(u,u,u) =F(u). (1.6) 

We assume that the zero equilibrium state always exists with .F(O) = and F'(0) = 1. This normal- 
ization allows us to represent the nonlinearity function in the form: 

f(u n -i,u n ,u n+ i) =u n + Q(u n -i,u n ,u n+ i), (1.7) 

where the linear part is uniquely normalized (parameter h can be chosen so that the linear term of 
u„_i + u n+ \ is cancelled) and the nonlinear part is represented by the function Q(u n -i,u n ,u n+ i). 
In addition, we assume that (i) F(u) and Q(u n -i, u n , u n+ i) are odd such that F(—u) = —F(u) and 
Q(—v, —u, —w) = —Q(v, u, w) and (ii) a pair of non-zero equilibrium points u + = —u_ ^ exists, 
such that 

F(u±) = 0, F'{u±) < 0, (1.8) 

while no other equilibrium points exist in the interval u G For instance, this assumption is 

verified for the 4 model il.31 with u± = ±1 and for the sine-Gordon equation (11.41 ) with u± = ±ir. 

We address the fundamental question of existence of travelling wave solutions in the discrete Klein- 
Gordon lattice (ll.lt . Since discrete equations have no translational and Lorentz invariance, unlike the 
continuous Klein-Gordon equation (11.2b . existence of travelling waves, pulsating travelling waves and 
travelling breathers represents a challenging problem of applied mathematics (see recent reviews in 
IS031 11.1051 1 

Our work deals with the travelling kinks between the non-zero equilibrium states u±. We are not 
interested in the travelling breathers and pulsating waves near the zero equilibrium state since the zero 
state is linearly unstable in the dynamics of the discrete Klein-Gordon lattice (ll.lt with F'(0) = 1 > 
0. Indeed, looking for solutions in the form u n (t) = e lKhn + Xn t we derive the dispersion relation for 
linear waves near the zero equilibrium state: 



It follows from the dispersion relation that there exists (h) > such that A 2 > for < k < (h). 
On the other hand, the non-zero equilibrium states u + and u- are neutrally stable in the dynamics of 
the discrete Klein-Gordon lattice dl.lt with F'(u±) < 0. We focus hence on bounded heteroclinic 
orbits which connect the stable non-zero equilibrium states u_ and u + in the form: 

Un(t) = 4>(z), z = hn — ct, (1.9) 

where the function <j>(z) solves the differential advance-delay equation: 

cV'(,) = + + *(*-*) + m + QWz _ hUizU(z + h)) . (i.io) 

We consider the following class of solutions of the differential advance-delay equation dl.lOt : (i) <j>{z) 
is twice continuously differentiable function on z G R; (ii) 4>{z) is monotonically increasing onzfl 
and (hi) 4>{z) satisfies boundary conditions: 

lim cj)(z)=u-, lim (f>(z) = u+. (1-H) 

It is easy to verify that the continuous Klein-Gordon equation (ll.2t with F(u) in (ll.8t yields a trav- 
elling kink solution in the form u = <fi(z), z = x — ct for |c| < 1. However, the travelling kink 
can be destroyed in the discrete Klein-Gordon lattice dl.lt . which results in violation of one or more 
conditions on <fi(z). For instance, the bounded twice continuously differentiable solution cj)(z) may 
develop non- vanishing oscillatory tails around the equilibrium states u± (TJ05 |. 

A recent progress on travelling kinks was reported for the discrete </> 4 model. Four particular spa- 
tial discretizations of the nonlinearity F(u) = u(l — u 2 ) were proposed with four independent and 
alternative methods IBT971 IS971 FZK99I IK03I . where the ultimate goal was to construct a family 
of translation-invariant stationary kinks for c = 0, that are given by continuous, monotonically in- 
creasing functions <j>[z) on z £ R with the boundary conditions dl.llt . Exceptional discretizations 
were generalized in [BOP05 DK Y05al IDKY05bll . where multi-parameter families of cubic polyno- 
mials f(u n -i, u n , w n +i) were obtained. It was observed in numerical simulations of the discrete </> 4 
model IIS97I that the effective translational invariance of stationary kinks implies reduction of radia- 
tion diverging from moving kinks. New effects such as self-acceleration were reported for some of the 
exceptional discretizations in I DKY05ai IDKY05B1 . Nevertheless, from a mathematically strict point 
of view, we shall ask if exceptional discretizations guarantee existence of true travelling kinks (hetero- 
clinic orbits) at least for small values of c. The question was answered negatively in [OPB05|, where 
numerical analysis of beyond-all-order expansions was developed. It was shown that bifurcations of 
travelling kink solutions from h = to small non-zero h do not generally occur in the discrete </> 4 
model with small values of c, even if the exceptional discretizations allow these bifurcations for c = 0. 
It was also discovered in [OPB05] that bifurcations of travelling kink solutions may occur for finitely 
many isolated values of c far from the limit c = 0. 
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The discrete sine-Gordon model was also subject of recent studies. Numerical computations were used 
to identify oscillatory tails of small amplitudes on the travelling kink solutions [EF90 SZE00|. These 
tails were explained with analysis of central manifold reductions (carried out without the normal form 
reductions) [FR05|. Exceptional discretizations of the nonlinearity F(u) = sin(u) were suggested 
by the topological bound method in [SW94] (see the review in [S99]) and by the "inverse" (direct 
substitution) method in [FZK99]. Simultaneously with the absence of single kinks in the sine-Gordon 
lattices, multiple kinks (between non-zero equilibrium points of 7r(mod27r)) were discovered with the 
formal reduction of the discrete lattice dl.lt to the fourth-order ODE problem in [ CKOO | and confirmed 
with numerical analysis of the differential advanced-delay equation in [ACR03]. 

Our work is motivated by the recent advances in studies of differential advanced-delay equations from 
the point of dynamical system methods such as central manifold reductions and normal forms (see 
pioneer works in [ 100 IK00 1). The same methods were recently applied in [ PR05 1 to travelling solitary 
waves in discrete nonlinear Schrodinger lattices near the maximum group velocity of linear waves. It 
was shown in [PR05] that non-existence of travelling solitary waves can be derived already in the 
truncated (polynomial) normal form. We shall exploit this idea to the class of travelling kinks in 
discrete Klein-Gordon lattices near the specific speed c = 1 and the continuous limit h = 0. This 
particular point conesponds to the quadruple zero eigenvalue on the central manifold of the dynamical 
system. A general reversible normal form for the quadruple zero eigenvalue was derived and studied 
in [195 ] (see also tutorials in the book [IA98 1). 

We note that the discrete Klein-Gordon lattice was considered in [IK00| but the nonlinearity F(u) was 
taken to be decreasing near u = 0, such that the quadruple zero eigenvalue was not observed in the list 
of possible bifurcations of travelling wave solutions (see Figure 1 in [IK00]). On the other hand, the 
quadruple zero eigenvalue occurred in the discrete Fermi-Pasta-Ulam lattice studied in [ 100 1 , where 
the symmetry with respect to the shift transformation was used to reduce the bifurcation problem to the 
three-dimensional center manifold. Since the reversible symmetry operator for the Fermi-Pasta-Ulum 
lattice is minus identity times the reversibility symmetry for the Klein-Gordon lattice, the reversible 
normal forms for quadruple zero eigenvalue are different in these two problems. We focus here only 
on the case of the discrete Klein-Gordon lattice. 

Our strategy is as follows. We shall develop a decomposition of solutions of the differential advance- 
delay equation into two parts: a four-dimensional projection to the subspace of the quadruple zero 
eigenvalue and an infinite-dimensional projection to the hyperbolic part of the problem. By using a 
suitable scaling, we shall truncate the resulting system of equations with a scalar fourth-order equation, 
which is similar to the one formally derived in [CK00|. We refer to this fourth-order equation as 
to the normal form for travelling kinks. We shall develop a numerical analysis of the fourth-order 
equation and show that no monotonic heteroclinic orbits from u_ to u + exist both in the discrete (f> A and 
sine-Gordon lattices. Rigorous persistence analysis of bounded solutions (heteroclinic orbits between 
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periodic solutions near u_ and iq_) is developed with the technique of center manifold reduction. Our 
main conclusion is that the differential advance-delay equation (ll.lOt has no monotonic travelling 
kinks near the point c = 1 and h = but it admits families of non-monotonic travelling kinks with 
oscillatory tails. 

It seems surprising that the truncated normal form is independent on the discretizations of the nonlin- 
earity function / (u n _ \ , u n , u n+ \ ) and the negative result extends to all exceptional discretizations con- 
structed in IISW94I 15331 IBT97I IFZK99I IKU51 IBOP05I IDKY05al IDKY05bl . Any one-parameter 
curves of the monotonic travelling kinks, which bifurcate on the plane (c, h) from the finite set of 
isolated points (c* , 0) (see numerical results in [OPB05 1), may only exist far from the point (1, 0). 

In addition, we shall explain bifurcations of multiple kinks from it_ = — it to u + = ir(2n — 1), 
n > 1 in the discrete sine-Gordon equation from the point of normal form analysis. These bifurcations 
may occur along an infinite set of curves on the plane (c, h) which all intersect the point (1,0) (see 
numerical results in [ ACR03 1). We shall also derive the truncated normal form from the discretizations 
of the inverse method reported in [FZK99 1 and show that it may admit special solutions for monotonic 
kinks, when the continuity condition (ll.6t is violated. 

We emphasize that our methods are very different from the computations of beyond-all-order expan- 
sions, exploited in the context of difference maps in [TTJ98 T00| and differential advanced-delay 
equations in IOPB 05 1. By working near the particular point c = 1 and h = 0, we avoid beyond-all- 
order expansions and derive non-existence results at the polynomial normal form. In an asymptotic 
limit of common validity, bifurcations of heteroclinic orbits in the truncated normal form can be stud- 
ied with beyond-all-order computations (see recent analysis and review in ITP05I ). 

Our paper has the following structure. Section 2 discusses eigenvalues of the linearization problem 
at the zero equilibrium point and gives a formal derivation of the main result, the scalar fourth-order 
equation. The rigorous derivation of the scalar fourth-order equation from decompositions of solutions, 
projection techniques and truncation is described in Section 3. Section 4 contains numerical analysis 
of the fourth-order equation, where we compute the split function for heteroclinic orbits. Persistence 
analysis of bounded solutions in the differential advance-delay equation is developed in Section 5. 
Section 6 concludes the paper with summary and open problems. A number of important but technical 
computations are reported in appendices to this paper. Appendix A contains the comparison of the 
fourth-order equation with the normal form from [I95|. Appendix B gives the proof of existence of 
center manifold in the full system, which supplements the analysis in [IK00|. Appendix C describes 
computations of the Stokes constant for heteroclinic orbits in the fourth-order differential equation 
from methods of [TP05 1. Appendix D discusses applications of the fourth-order equation to the inverse 
method from [FZK99|. 
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2 Resonances in dispersion relations and the scalar normal form 



When the differential advance-delay equation dl.lOt is linearized near the zero equilibrium point, we 
look for solutions in the form <j){z) = e Xz , where A belongs to the set of eigenvalues. All eigenvalues 
of the linearized problems are obtained from roots of the dispersion relation: 

D(A;c,h) = 2(coshA-l) + /i 2 -c 2 A 2 = 0, (2.1) 

where A = Ah is the eigenvalue in zoomed variable C = f ■ We are interested to know how many 
eigenvalues occur on the imaginary axis and whether the imaginary axis is isolated from the set of 
complex eigenvalues (that is complex eigenvalues do not accumulate at the imaginary axis). Imaginary 
eigenvalues A = 2iK of the dispersion relation (I2.lt satisfy the transcendental equation: 

h 2 

sin 2 K = — + c 2 K 2 . (2.2) 

Particular results on roots of the transcendental equation (12.21 ) are easily deduced from analysis of 
the function sin(x) in complex domain. Due to the obvious symmetry, we shall only consider the 
non-negative values of c and h. 

• When c = and < h < 2, equation ( 12.21 ) has only simple real roots K, such that all eigenval- 
ues A are purely imaginary and simple. All real roots K become double at h = and h = 2. 

• When h = and 1, a double zero root of K (a double zero eigenvalue A) exists. When 
< c < 1, finitely many purely imaginary eigenvalues A exist (e.g. only one pair of roots 
K = ±P exists for cq < c < 1, where cq rj 0.22), while infinitely many roots are complex 
and distant away the imaginary axis. When c > 1, all non-zero eigenvalues A are complex and 
distant away the imaginary axis. 

In the general case of h ^ and c / 0, the transcendental equation d2.2t is more complicated but it 
can be analyzed similarly to the dispersion relation considered in [100, IK00]. 

Lemma 2.1 There exists a curve h = /i*(c), which intersects the points (1, 0) and (0, 2) on the plane 
(c, h), such that for < c < 1 and < h < h*{c) finitely many eigenvalues A of the dispersion 
relation i2.lt are located on the imaginary axis and all other eigenvalues are in a complex plane distant 
from the imaginary axis. The curve h = /i*(c) corresponds to the 1:1 resonance Hopf bifurcation, 
where the set of imaginary eigenvalues includes only one pair of double eigenvalues. 

Remark 2.1 The 1:1 resonant Hopf bifurcation is illustrated on Figure 1. Figure 1(a) shows a graph- 
ical solution of the transcendental equation \2.2\ for h = h*(c). Figure 1(b) shows the bifurcation 
curve h = /^(c) on the plane (c, h). The curve h = for < c < 1 corresponds to the double zero 
eigenvalue in resonance with pairs of simple purely imaginary eigenvalues. 



6 



Proof. Let A = p + iq and rewrite the dispersion relation in the equivalent form: 

c 2 (q 2 -p 2 ) + h 2 + 2(coshpcosg - 1) = 
c 2 pq — sinhpsing = 0. 

It follows from the system with c / that the imaginary parts of the eigenvalues is bounded by the 
real parts of the eigenvalues: 

q 2 <p 2 + 4 cosh2 ?- ( 2 -3) 
c z 2 

Therefore, if complex eigenvalues accumulate to the imaginary axis, such that there exists a sequence 
(PmQn) with lim n -+ 00 p n = and lim n ^oo q n = q*, the accumulation point (0,g*) is bounded. 
However, since D(A; c, h) is analytic in A, the dispersion relation (I2.lt may have finitely many roots 
of finite multiplicities in a bounded domain of the complex plane. Therefore, the accumulation point 
(0, q*) does not exist and complex eigenvalues are distant from the imaginary axis. By the same reason, 
there exist finitely many eigenvalues on the imaginary axis A € iR. 

Looking at the double roots at K = ±P, we find a parametrization of the curve h = h* (c) in the form: 

c 2 = smP ^ osP 7 h 2 =4smP( S mP-PcosP). (2.4) 

A simple graphical analysis of the transcendental equation (12.2b shows that the double roots at K = 
±P are unique for < P < ^, when the intersection of the parabola and the trigonometric function 
occurs at the first fundamental period of sin 2 -function. At P = (at the point (c,h) = (1, 0)), the pair 
of double imaginary eigenvalues A = ±2iP merge at A = and form a quadruple zero eigenvalue. 
At P = 5 (at the point (c, h) = (0, 2)), a sequence of infinitely many double imaginary eigenvalues 
exists at A = m{l + 2n), n G Z. ■ 



We will be interested in the reduction of the differential advanced-delay equation dl.lOt at the partic- 
ular point (c, h) = (1, 0). Let e be a small parameter that defines a point on the plane (c, h), which is 
locally close to the bifurcation point (1,0): 

c 2 = 1 + e 7 , h 2 = e 2 r. (2.5) 

Let A = v^Ai be the scaled eigenvalue near A = 0, so that the Taylor series expansion of the disper- 
sion relation (I2.lt allows for a non-trivial balance of e-terms. Truncating the Taylor series beyond the 
leading order 0(e 2 ), we obtain the bi-quadratic equation for the scaled eigenvalue Ai: 

^Aj - 7 Aj + t = 0. (2.6) 

The 1 : 1 resonance Hopf bifurcation corresponds to the point, where the roots of the bi-quadratic equa- 
tion (I2.6t are double and purely imaginary. This point occurs when r = 3^f 2 , which agrees with the 
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Figure 1 : Left: The graphical solution of the transcendental equation d2.2t for P = ?, where (c, h) are 
given by A2.4I ). Right: The bifurcation curve for 1:1 resonant Hopf bifurcation h = h*(c) and double 
zero resonance bifurcation h = 0. 

Taylor series expansion of the system (12.41) in powers of P at the leading orders of P 2 and P 4 . As a 
result, the bifurcation curve h = h* (c) has the leading order behavior: 

K{c) = \/3(l-c 2 )+o(l-c 2 ). (2.7) 

The normal form for travelling kinks, which is the main result of our paper, can be recovered with a 
formal asymptotic expansion of the nonlinear differential advance-delay equation dl.lOt . Let C = f 
and rewrite the main equation ( 11.101 ) in the form: 

cV'CO = <KC + 1) - 2^(0 + <XC - 1) + h 2 4>{Q + h 2 Q(cf>(( - 1), 0(C), HC + !))• (2-8) 

Let (c, h) be the perturbed bifurcation point in the form d2.5t . Let 0(C) be a smooth function of a 
slowly varying variable Ci = \/eC- Expanding 0(Ci ± y^e) in the Taylor series in y/e and truncating 
by terms beyond the leading order e 2 , we obtain the scalar fourth-order differential equation: 

J_0(iv) (Cl) _ 7 0"(d) +rF(0(Ci)) = 0, (2.9) 

where we have used the representations ( 11.61 ) and dl.71 ). The linearization of the nonlinear ODE 
( I2.9t near the zero equilibrium point recovers the bi-quadratic dispersion relation (I2.6I) . The non- 
linear equation ( 12.91 ) has the equilibrium points U-, 0, and u + , which are inherited from the equilib- 
rium points of the Klein-Gordon equation (11.21 . Our main examples will include the 4 model with 
F(u) = u(l — u 2 ), u± = ±1 and the sine-Gordon equation with F{u) = sin(u), u± = ±7r. 
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3 Decompositions, projections and truncation 



We shall derive the normal form equation (I2.9t with rigorous analysis when the solution of the dif- 
ferential advance-delay equation (I2.8t is decomposed near a quadruple zero eigenvalue. We adopt 
notations of [IKOQl (see review in IIJ05I ) and rewrite the differential advance-delay equation as an 
infinite-dimensional evolution problem. We shall work with the scaled (inner) variable C = § > where 
the differential advance-delay equation dl.lOt takes the form d2.8t . Let p be a new independent vari- 
able, such that p £ [-1, 1] and define the vector U = (E7i(0, U 2 ((), U 3 ((,p)) T , such that 



U 



0(C), u 2 = <p'(0, u 3 = 0(C + P ). 



(3.1) 



It is clear that Us((,p) = Ui(( + p) and U^{C,, 0) = U\{C). The difference operators are then defined 
as 



U 3 ((, ±1) = 8±U 3 (C,p) = 0(C ± 1) = 6±Ux(0 = C/a(C ± 1). 
Let V and 7i be the following Banach spaces for U = (U\, U 2 , U 3 {p)) T , 

v = {(u l ,u 2 )e^ 2 ,u^c l ([-iM^),um = Ui}^ 
n = {(u u u 2 ) e« 2 , c/ 3 gc ([-i,i],r)}, 



(3.2) 

(3.3) 
(3.4) 



with the usual supremum norm. We look for a smooth mapping £ ^ U(0 in C°(W;T>), which 
represents classical solutions of the infinite-dimensional evolution problem: 



A U = £ U + ^ Mo (U), 
where C Cj h and Mq(U) are found from the differential advance-delay equation (12.81) : 



(3.5) 



C 



c,h 



\ 



X(5+ + 5- 







On 



and 



M (U) = {Q,Q{5-U z ,U 1) 5 + Uz)^) T . 



(3.6) 



(3.7) 



The linear operator C C: h maps V into TL continuously and it has a compact resolvent in TL. The 
nonlinearity Mo(U) is analytic in an open neighborhood of U = 0, maps TL into V continuously, and 
||M (U)b = O(||U|^). 

The spectrum of C c ^ consists of an infinite set of isolated eigenvalues of finite multiplicities. By 
virtue of the Laplace transform, eigenvalues of the linear operator C Cj h are found with the solution 
U(C,p) = (1, A, e Ap ) T e A ^, when the linear problem U'(£) = C C ^XJ is reduced to the dispersion 
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relation (I2.lt . that is D(A; c, h) = 0. We are particularly interested in the bifurcation point c = 1 and 
h = 0, when 

£>i(A) = D(A;1,0) = 2(cosh A - 1) - A 2 . (3.8) 

The transcendental equation D\ (A) = has the quadruple zero root and no other root in the neigh- 
borhood of A G iR. The four generalized eigenvectors of the Jordan chain for the zero eigenvalue, 
CifiXJj = Uj-i, j = 1, 2, 3, 4 with U_i = 0, are found exactly as: 





/ 


1 






( o\ 




f 


o ^ 




/ 









U = 









Ui = 


1 


u 2 = 







u 3 = 









(3.9) 




V 


1 


J 




[p ) 




\ 


\v 2 ) 




\ 


6^ 







The dynamical system (13.51) has the reversibility symmetry S, such that 



where 



su 



£ c , h 5U + ^M (5U), 



(3.10) 



(3.11) 



(u 1 (C),-u 2 (C),u 3 (c,-p)) t 

and the symmetry property (ll.5t has been used. Since both the linear operator d3-6l > and the nonlinearity 
i3.lt anti-commute with the reversibility operator (13. Ill) , the standard property of reversible systems 
holds: if U(£) is a solution of d3.5t for forward time £ > 0, then <SU(— £) is a solution of d3.5t for 
backward time ( < (see [LR98 1 for review). Applying the reversibility operator S to the eigenvectors 
d3-9l > at the bifurcation of the quadruple zero eigenvalue, we observe that 



SVj = (-lyVj, j = 0,1,2,3. 



(3.12) 



This bifurcation case fits to the analysis of [195 1, where the reversible normal form was derived for the 
quadruple zero eigenvalue. Appendix A shows that the general reversible normal form from [193 is 
reduced to the normal form equation i2.9i by appropriate scaling. However, we notice that this result 
cannot be applied directly, since it is only valid in a neighborhood of the origin, while in the present 
case, the solutions we are interested in are of order of 0(1). 

In order to study bifurcation of quadruple zero eigenvalue in the reversible system (I3.5I) . we shall define 
the perturbed point (c, h) near the bifurcation point (1, 0) with an explicit small parameter e. Contrary 
to the definition < I2.5I >. it will be easier to work with the parameters (7, r), defined from the relations: 



1 



1-67, 



,-2 



S 2 T, 



(3.13) 



where e is different from e used in (I2.5t . The dynamical system (13.51) with the parametrization (I3.13t 
is rewritten in the explicit form: 



d 



U = £1 U + e 7 Li(U) + £ 2 tL 2 (U) + £ 2 rM (U), 



(3.14) 
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where Cifl follows from (I3.6t with (c, h) = (1, 0), Mo(U) is the same as (I3.7I) . and the linear terms 
Li ;2 (U) are 

Li(U) = (0,2U 1 -(5 + + 5-)U 3 ,0) T , L 2 (U) = (0,C/i,0) T . (3.15) 

Bounded solutions of the ill-posed initial-value problem ( 13.141 ) on the entire axis ( S tare subject of 
our interest, with the particular emphasis on the kink solutions (heteroclinic orbits between non-zero 
equilibrium points u_ and u + ). These bounded solutions can be constructed with the decomposition 



of the solution of the infinite-dimensional system d3.14t . projection to the finite-dimensional subspace 
of zero eigenvalue and an infinite-dimensional subspace of non-zero eigenvalues, and truncation of 
the resulting system of equations into the fourth-order differential equation d2.9t . Since the original 
system is reversible, the fourth-order equation must inherit the reversibility property. 

The techniques of decompositions, projections and truncation rely on the solution of the resolvent 
equation (AT - £i,o)U = F, where U G V, F G H and A G C. When A is different from the roots 
of Di(A), the explicit solution of the resolvent equation is obtained in the form: 

C/x = -[DiCA)]- 1 ^), 

U 2 = At/! -F u 

U 3 (p) = £/ie Ap - [ P F 3 (s)e A ^- s Us, (3.16) 
J o 

where 

F(A) = AFi + F 2 - J (F 3 (s)e A ^ - F 3 (- S )e~^ 1 - S ^ ds. (3.17) 

The quadruple zero eigenvalue appears as the quadruple pole in the solution TJ(p; A) of the resolvent 
equation near A = 0. Let us decompose the solution of the dynamical system (13.141) into two parts: 

U = X + Y, (3.18) 

where X is a projection to the fourth-dimensional subspace of the quadruple zero eigenvalue, 

X = A(C)U + B(C)Ui(p) + C(C)U 2 (p) + D(()XJ 3 (p) (3.19) 

and Y is the projection on the complementary invariant subspace of >Ci o- We notice in particular that 

EMC) = M0 + Yi(Q (3.20) 

U 3 ((,p) = A(()+ P B(() + \p 2 C(0 + lp 3 D(() + Y 3 ((,p). (3.21) 

2 b 



With the decomposition d3.18t - d3.19l ). the problem d3.14t is rewritten in the form: 

^-Y-£ lt0 Y = F(A,B,C,D,Y), (3.22) 
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where 



Fi 
F 2 
F 3 



B-A! 

C-B'-£ 7 (C- 2Yt + (S + + 5-)Y 3 ) + e 2 t{A + Q + Fi 

B — A! + p(C — B') + \p\D - C) + \p\-D'), 

2 o 



and 



q = q[a-b + \c-\d + S~Y 3 , A + Yi, A + B + ~C + \d + 5+y 3 

2 b 2d 



We notice that p is a dumb variable in (I3.22I) . Since the vector Y is projected to the subspace of non- 



zero eigenvalues of £1.0, the right-hand-side of the "inhomogeneous" problem (I3.22t must be orthog- 
onal to eigenvectors of the adjoint operator £* associated to the zero eigenvalue. Equivalently, one 
can consider the solution of the resolvent equation (I3.17t and remove the pole singularities from the 
function [Di(A)] _1 F(A) near A = 0. See p. 853 in [ 100 1 for the projections on the four-dimensional 
subspace and, in particular, the formula: 



/ \ 



-1 







\ 



= F (p), 



(3.23) 



\ J \ - hp 2 ) ) 

where Vh is a projection operator to the space for Y. This procedure results in the fourth-order 
differential system for components (A, B,C,D) 6R 4 in £: 



A' 
B' 
C 
D' 



B 
C 
D 



~£ 7 (C - 2Yt + (6+ + r )Y 3 ) + ~ £ 2 r (A + Q + Y x ) 



12ej (C - 2Yi + (5+ + <T)Y 3 ) - 12e 2 r (A + Q + Yi). 



(3.24) 



While the system (I3.22t has the reversibility symmetry S in ( 13.1 It . the system (13.241) has the reduced 
reversibility symmetry, 



S A = (A,-B,C, -D) T , 



(3.25) 



where A = (A, B,C, D) T . The closed system of equations ( 13.221 ) and ( 13.241 ) can be rewritten by 
using the scaling transformation: 

A = MCi), B = Vi<MCi), C = eMCi), D = e 3 / 2 Md), Y = e 3 / 2 ^((,p), (3.26) 
where £i = yfeC,. The four-dimensional vector cf> G M 4 satisfies the vector system: 



03 



/ 1 
0010 
0001 
\ 



03 

V ^ 4 / 



+ 



/ \ 



12 



(3.27) 
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where g is a scalar function given by 

g = 7 (0 3 - 2yfhh + + <nV> 3 ) - t($i + Q + e 3/2 ^i) 

and 

Q = qUi- Vlfo + ^efa - ^ 3/2 04 + V^S'ips, <j>\ + e 3/2 ipi, 

01 + ^<p2 + ^£03 + £^04 + ^S+lfaj ■ 

The infinite-dimensional vector if) G P satisfies the "inhomogeneous" problem: 

A.oV = ^F (p), (3.28) 



where Fo is given in d3.23t . The truncation of the fourth-order system (13.271) at e = recovers 
the scalar fourth-order equation d2.9t . The set of equilibrium points of the coupled system (13.27 1 >— 
(I3.28t contains the set of equilibrium points of the fourth-order equation d2.9t . since the constraints 
02 = 03 = 04 = and = reduce the coupled system to the algebraic equation <fii + Q (0i , 0i , 0i ) . 
that is F(4>i) = 0. Our central result is the proof of existence of the center manifold in the coupled 



system (I3.27t - d3.28l l. 

Theorem 1 Fix M > and let e > be small enough. Then for any given 

IMI<M, (3.29) 
there exists a unique solution ip G C°(R, V) of the system \3.28\ . such that 

< sfeK, (3.30) 



where the constant K > Ois independent of e. Moreover, when (f)((i) is anti-reversible with respect to 
So, i.e. 5o0(— Ci) = — <^(Ci)> men V'(C) i s anti-reversible with respect to S, i.e. Sij}{—Q = —■0(C). 
If the solution 4>(Ci) tends towards a T -periodic orbit as Ci — * =too, then the solution tends 



towards a T j ^fit-periodic orbit of the system \3. 28\ as Q — > ±oo. 

In this theorem, C° (R, 6) denotes the Banach space (sup norm) of continuous and bounded functions 
on R taking values in the Banach space £. 

The proof of existence of the center manifold for the class of reversible dynamical system (I3.5t is 
developed in Lemmas 2-4 of HIK00I with the Green function technique. Applications of the center 
manifold reduction in a similar content can be found in lld0l llJ()5i rFR05 1. We note that the case of the 
quadruple zero eigenvalue was not studied in any of the previous publications (e.g. the pioneer paper 
[IK00| dealt with the case of two pairs of double and simple imaginary eigenvalues). We develop 



an alternative proof of existence of center manifold in the coupled system d3.27t - d3.28t with a more 
explicit and elementary method in Appendix B. 



13 



4 Numerical analysis of the truncated normal form 



Eigenvalues of the linear part of the fourth-order equation A2.91 . that are roots of the bi-quadratic 
equation J2.6t . are shown on Figure 2. Due to the definition (I3.13l l. it makes sense to consider only the 
upper half plane of the parameter plane (7, r). The curve Ti = {7 < 0, r = 0} corresponds to the 
double zero eigenvalue that exists in resonance with a pair of purely imaginary eigenvalues. The curve 
T2 = {7 < 0, r = 37 2 } corresponds to the reversible 1:1 resonance Hopf bifurcation h = /i*(c) in the 
limit i2.lt . These two curves confine the parameter domain of our interest, where a heteroclinic orbit 
between two non-zero equilibrium points u_ and u + of the nonlinearity function F(u) may undertake 
a bifurcation. 



X 




(0,0) 



Figure 2: Eigenvalues of the scalar normal form ( 12.91 ) on the parameter plane (7, r). The curves ^,1,2,3 
mark various bifurcations of co-dimension one. 

A complete classification of various bifurcations of co-dimension one in the normal form i2.9i can 
be found in [195]. Homoclinic orbits may bifurcate above the curve T2 and may exist in the domain 
D = {t > 37 2 , 7 < 0} U {r > 0, 7 > 0}. Homoclinic orbits correspond to localized solutions to the 
zero equilibrium state and they are not studied in this paper. We focus hence on bounded heteroclinic 
orbits which connect the non-zero equilibrium states u_ and u + . 

Let us rewrite the scalar fourth-order equation d2.9t in the normalized form, 

0(*O + a( j)» + F(<t)) = 0, (4.1) 
where 4> = 4>{t), t = (12r) 1//4 (j, and a = — ^| 7 - The domain between the curves T2 and Ti 
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corresponds to the semi-infinite interval a > 2. We shall consider separately the discrete <^> 4 model 
with F(4>) = 0(1 — (f) 2 ) and the sine-Gordon model with F(4>) = sin(0). 

When F{4>) = 4>(1 — <ft 2 ), there are two non-zero equilibrium points u± = dbl. Linearization of the 
scalar equation (I4.lt near non-zero equilibrium points gives a pair of real eigenvalues (Ao, — Ao) and a 
pair of purely imaginary eigenvalues (iu>o, —iujq) for any a G R. Our numerical algorithm is based on 
the consideration of the unstable solution 4> u (t), such that 

lim (j> u (t) = -l, lim ((P u (t) + l)e- Xot = c , (4.2) 

t— >— oo t^— oo 

where co is arbitrary positive constant, which is related to the translation of the unstable solution. Let 
to(co) be a zero of the unstable solution <f) u (t), if it exists. Since the ODE problem (I4.lt with odd 
function F (—(/>) = —F((/>) has the reflection symmetry <fi h-> —(f), we look for a bounded heteroclinic 
orbit as an odd function of t. Therefore, the distance between stable and unstable solutions can be 
measured by the split function K = 0"(to(co)), such that the roots of K give odd heteroclinic orbits. 
Although the split function K is computed at to, which depends on cq, it follows from the translational 
invariance of the problem (I4.lt that K is independent of cq. (It may depend on the external parameter 
a of the problem (I4.lt .) The first zero to corresponds to the monotonic heteroclinic orbit when K = 
0, while subsequent zeros correspond to non-monotonic heteroclinic orbits. We shall here consider 
monotonic heteroclinic orbits. 

The unstable solution of the scalar equation (I4.lt is numerically approximated by the solution of the 
initial-value problem with the initial data: 





m 


\ 
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-1 
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i 




^(0) 











+ c 




Ao 


















(/>"(0) 













Ao' 
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0"'(O) 


/ 


V 





) 




V 


Aq 



where Ao > and < cq <C 1. The solution of the initial- value problem is stopped at the first 
value to > 0, where </>(to) = and the split function K is approximated as K = (//'(to). Besides the 
external parameter a, the numerical approximation of the split function K may depend on the shooting 
parameter Co and the discretization parameter dt of the ODE solver, such that K = K(a; co, dt). 

The graph of K versus a for a fixed set of values of co and dt is shown on Figure 3(a). It is clearly seen 
that the split function K(a) is non-zero for finite values of a. It becomes exponentially small in cr _1 
as a — > oo. In the beyond-all-order asymptotic limit, computations of Stokes constants give equivalent 
information to the computation of the split function K(a) (see BTTJ98I T00|). We apply this technique 
in Appendix C to show that the Stokes constant for the scalar fourth-order equation (I4.lt is non-zero 
in the limit a — > oo. 

Behavior of the split function K for a fixed value of a is analyzed on Figures 4(a,b) as the values of 
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dt and cq are reduced. When dt is sufficiently reduced, the value of K converges to a constant value 
monotonically (see Figure 4(a)). When Co changes, the value of K oscillates near a constant level. 
The amplitude of oscillations depend on the value of dt and it becomes smaller when the value of dt 
is reduced (see Figure 4(b)). Oscillations on Figure 4(b) can be used as the tool to measure the error 
of numerical approximations for the split function K{a). The maximum relative error is computed 
from the ratio of the standard deviation of the values of K(a) to the mean value when the parameter cq 
varies. The graph of the relative error is plotted on Figure 3(b) versus parameter a. The relative error 
increases for larger values of a since K(a) becomes exponentially small in the limit a — > oo. 




Figure 3: Left: The graph of the split function K(a) for Co = 0.00001 and dt = 0.005 for kink 
solutions of the ODE (14.11) with F(<p) = (p(l — (j) 2 ). Right: relative error of numerical approximations 
of K(a) for dt = 0.005 under variations of cq. 




Figure 4: Left: Convergence of the split function K(a) as dt — > for a = 5 and Co = 0.00001. Right: 
Oscillations of the split function K{a) as cq —> for a = 5 and dt = 0.005. 
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When F{4>) = sin(0), there are two non-zero equilibrium points u± = ±7T. Linearization of the 
scalar equation (I4.lt near non-zero equilibrium points still has a pair of real eigenvalues (Aq, — Ao) 
and a pair of purely imaginary eigenvalues (iuo, —iu>o) for any a G E. Therefore, we adjust the same 
algorithm to the sine-Gordon model. Figures 5(a) displays the graph of the split function K versus a 
for a fixed set of values of cq and dt, which Figure 5(b) shows the maximum relative error of numerical 
approximations. The split function K(a) displays the same behavior as that for the 4 model. 




Figure 5: Left: The graph of the split function K(a) for Co = 0.00001, dt = 0.005, and different 
values of a for kink solutions of the ODE (I4.lt with F{4>) = sin(</>). Right: relative error of numerical 
approximations of K(a) for dt = 0.005 under variations of cq. 

Since the discrete sine-Gordon lattice has multiple equilibrium states at <fi = ttu, n € Z, one can con- 
sider multiple monotonic kinks between several consequent equilibrium states [CK00|. In particular, 
we shall consider the double kink, which represents a bound state between the heteroclinic orbit from 
u_ = —it to u_|_ = 7T and the heteroclinic orbit from u + = it to u + = 3tt. 

Double kinks were considered recently in a fourth-order differential model in [CK00| and in the differ- 
ential advance-delay equation in [ACR03|. We note that the fourth-order model in [CK00| is obtained 
in the continuous limit of the discrete Klein-Gordon lattice (ll.lt . such that the fourth-order derivative 
term is small compared to the second-order derivative term. However, numerical results of [CK00| 
were derived in the case when both derivatives are comparable. Asymptotic and direct numerical re- 
sults showed that the families of double kinks intersect the bifurcation point (c,h) = (1,0), where our 
normal form (I2.9t is applicable. Thus, we give a rigorous explanation of existence of double kinks in 
the discrete sine-Gordon lattice from the theory of normal forms. 

We modify the numerical algorithm by considering the initial- value problem for the normalized fourth- 
order equation d4.lt with the initial values (I4.3t . Let to be defined now from the equation <j)(to) = u + , 
and the split function K(a) be defined by K = <j)"(to). When K = 0, the bounded solution <p u (t) 
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defines a double kink solution with the properties: lim^_oo <p u (t) = u_, lim^oo <p u (t) = u + and 
4> u (t — to) = —4> u {to — t) + 2tt. The graph of K versus a on Figure 6 shows infinitely many zeros 
of K(a), which correspond to infinitely many double kink solutions (only four such solutions were 
reported in JCKOOI). 



4, 

K(0) 

3 



2 







-1 



-2 




Figure 6: The graph of the split function K(a) for double kink solutions of the ODE (I4.ll) with 

F(4>) = sin(0). 

In addition, we mention that the numerical method of computations of the split function can be tested 
for other heteroclinic orbits, such as triple monotonic kinks between (— tt, tt), (7r,37r), (37t,57t) and 
triple non-monotonic kinks between (— 7T, tt), (tt, —it), (— 7r, it). It can also be used for construction 
of homoclinic orbits (see ITP05I ). The accuracy of computations of the split function is only limited 
to the accuracy of the ODE solver for the scalar fourth-order differential equation ( 14. It . 

5 Persistence analysis of bounded solutions 

We consider the truncation of the system (I3.27t in the form: 

= N(<£), (5.1) 
where <fi = (<^, 4>2, 4>3, 4>4) T and the vector field N : M 4 i— > M 4 is given explicitly in the form 

N(0) = (0 2 , 03, 04, 12 7 03 - 12tF(0!)) t . (5.2) 
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Let us use t = Ci in notations of this section. The truncated system (I5.lt exhibits a number of useful 
properties: 



PI Since F(fa) is odd in fa £ R, the vector field N(0) is odd in G R 4 

P2 The system is conservative, such that div N = 0. 

P3 The system is reversible under the reversibility symmetry So in d3.25t . 

P4 The set of fixed points of ±So is two-dimensional. 

P5 The system has two symmetric equilibrium points u± = (u±, 0, 0, 0) T . 

P6 Since F'(u±) < 0, the linearized operator Z)N(u±) has two simple eigenvalues on the imagi- 
nary axis and two simple real eigenvalues for any 7 £ R and 

We shall add a generic assumption on existence of a heteroclinic connection between periodic solutions 
near the non-zero equilibrium points. As it follows from numerical results in Section 4, we can not 
assume the existence of a true heteroclinic orbit between the two equilibrium points u± since the true 
heteroclinic orbits are not supported by the nonlinear functions of the discrete fa and sine-Gordon 
equations. However, a numerical shooting method implies an existence of a one-parameter family of 
anti-reversible heteroclinic connections between periodic solutions. Indeed, an initial-value problem 
for an anti-reversible solutions fat) of the truncated system ( 15. II ) (such that S^fa—t) = — fat)) has 
two parameters, while there is only one constraint on the asymptotic behavior of fat) as t — > 00. 
An orthogonality constraint must be added to remove an eigendirection towards a one-dimensional 
unstable manifold associated to the real positive eigenvalue Ao in the linearization of u + . It is expected 
that a one-parameter family of solutions of the initial-value problem satisfies the given constraint. We 
will make a generic assumption on existence of the one-parameter family of anti-reversible solutions 
of the truncated system (15.11) and we intend to prove that such solutions persist for the full system 



Assumption 1 For a fixed set of parameter (7, r) and a fixed nonlinear function F{fa), there exists 
a smooth family of solutions (f> a (t), such that 



0« = N(0 a ), 



(5.3) 



and 



•So^a(-t) 



4> a {t)- 



(5.4) 
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When t — » ±00, the family 4> a (t) tends towards u~(t), where u~(t) are periodic solutions of the 
system f l5.il ) near the equilibrium points u±, such that 

5 u-(-i) = -u+(i) (5.5) 

and 

S u~(-t- S a ) = u~(t- S Q ). (5.6) 
We assume that T a (period at infinity) and 5 a (phase shift) are smooth functions of parameter a. 

Let us choose a particular heteroclinic connection <p Q . We define a linear operator C associated with 
the linearized system d5.lt near 4> (t): 

C(t)cp(t) = cp(t) - DN(<t> (t))<p(t). (5.7) 

If the periodic solutions (t) have a sufficiently small radius, the spectrum of the linearized system 
d5.7b on the periodic solutions Uq (t) is predicted from that at the equilibrium points u^. Our aim is 
to study the invertibility of C in the space of continuous functions on t € M, which are asymptotically 
periodic as t — > ±cxd. In what follows, we prove the following result, under an additional assumption 
that the spectrum of C is generic (free of bifurcation). The Fredholm Alternative Theorem would 
need to be used in the case when the assumption fails and the kernel of C includes an anti-reversible 
decaying solution. 

Assumption 2 There is no anti-reversible solution of the homogeneous linearized system that decays 
to zero at both infinities, i.e. such that C(t)ip(t) = 0, So<p(—t) = —<p(t), and limt_ >r i- 00 <p(t) = 0. 

Theorem 2 Suppose As sumptions [I]and\2\hold. Let F £ C?(R, M 4 ) be a function that is reversible 
with respect to Sq and tends towards a To-periodic function as t — > ±00. Then, there exists a unique 
solution (p € C^(1R,]R 4 ) of the system 

C(t)<p(t) = F(t), (5.8) 

which is anti-reversible with respect to Sq and tends towards a To-periodic function as t — > ±00. The 
linear map F 1— ► cp is continuous in CP(K, M 4 ). 

Proof. We shall first consider the fundamental matrix for the time-dependent linear system, 

cp = DN(0o(*))¥>- (5-9) 
Due to properties of the truncated system d5.lt . we have the following relationship, 

cS OJ DN(0 o (t))y> = -DN(cf> (-t))S cp. (5.10) 
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By differentiating (15.31) with respect to t, we obtain the first eigenvector of C{t): 

£(t)0 o (i)=O. (5.11) 

By differentiation (I5.3t with respect to a, we can obtain another eigenvector of C(t). However, since 
the period depends on a, the second eigenvector is non-periodic at infinity but contains a linearly 
growing term in t. So, let us define a new "time" r as 

t = k a T, k a = T a /T , ko = 1. (5.12) 

Then 4> a (r) = 4> a (t) has the period To as r — > ±oo and satisfies 

-^-0 Q = fe a N(0J. (5.13) 

Differentiating this identity with respect to a and denoting </> = <9 a <^ a | Q =o we obtain the equation, 

J^o = DN(0o)$o + *iN(0 o ), ( 5 - 14 ^ 

where 

fcx = 5 a A;| a= o. (5.15) 
This explicit computation shows that O (*) * s tne generalized eigenvector of the operator C(t), 

C(t)fa(t) = hio®, (5.16) 

which is bounded and To-periodic at infinity as the first eigenvector 4>o(t). Moreover, we note the 
symmetry properties: 

SoM-t) = o (t), (5.17) 

S 4> (-t) = -0 o (t). (5.18) 



The two eigenvectors o (i) and <fi (t) correspond to the "neutral" directions associated with the fami- 

Iq and Uq 



lies of periodic solutions u„ (t) near the equilibria u±. Indeed, similarly we denote by Uq and Uq the 



To— periodic solutions of the limiting system with periodic coefficients: 

These two modes are the eigenvector and generalized eigenvector belonging to the double zero Flo- 
quet exponent associated with each periodic solution Uq (i). Since these periodic solutions are close 
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enough to the equilibria u± , and thanks to perturbation theory, we know that there exist two real Flo- 
quet exponents (Aq , — Aq ) for each solution Uq . Let us denote by C±(*) and X±(t) the To-periodic 
eigenvectors (defined up to a factor) such that 



: C + + A C+ = £N(u 5 )^ 



d_ 

|C--A C- = £»N(u 5 )C- 

d_ 

It 



X--KX- = DN(u+)x_. 



We justify below that Aq^ = A . Since the function SoC+( — t) satisfies 

J t S oC + (~t) - A SoC+H) = DN(Sou 5 (-t))SoC + (-t) (5-19) 
and because of (15.61) we can choose £_ such that 

S'oC + (-*-<5o) = C-(t-*o). (5-20) 
Now, since we have ( 15.51) and since DN(<fi (t)) is even in 4> (t), we deduce that 

£>N(u+(i)) = £>N(-S (uq (-*)) = L»N(uo (t - 26a)), 
such that Aq = Aq~. Finally we can choose the vector functions £± and x± sucn that 

X+it + h) = C+(*-^o) 

X-(t + 6 ) = C-(t-6 ) (5.21) 
SoX+i-t + h) = X-(t + 6 ). 

Since C±(t)e ±Xot and x±(t)e ±Xot are solutions of the limiting linear system, we introduce the unique 
solutions cf> 1 2 {t) and tp 1 2 (t) of the system i5.9i on t E M, such that 

<^iW ~ C + (i)e Ao ' as f -> -oo 
fait) ~ C-(*)e _Ao * as t^-cx) 

~ X+(*)e Aot as t -» oo 
?/>i(f) ~ X_0)e~ A °* as i -> oo. 

Now we have 

j^oM-t) = DN(SoM-t))SoM-t) = DN(Mt))SoM~t) 

and 

S <f>i(-t) ~ S'oC + (-i)e~ Ao * as t -» oo, 
22 



where 

S C+(-t)e- Xot = C_(t - 2S )e- Xot = x_(*)e~ A °'- 

Hence 

5 o 0i(-t) = ViO), (5-22) 

and, in the same way 

S o 2 (-t) = 2 (*)- ( 5 - 23 ) 

The two pairs of linearly independent solutions ((j> 1 , <p 2 ) and (tpi, ij} 2 ) cannot contain components in 
4> Q and 4> thanks to the behavior at infinity, hence there exist constants (a, b, c, d) such that 

V>! = a<t> x + bc/> 2 

iP 2 = c0 x + d<f> 2 . (5.24) 





Iterating once the relations (15.22b and (15.23b . we find that 

/ a 2 + bc b(a + d) \ 
\ c(a + cf) d 2 + bc J 

If either 6 or c is non-zero, then d = —a and a 2 + be = 1. If 6 = c = 0, then a 2 = d 2 = 1. The choice 
a = — d = ±1 is included in the previous case d = —a. The choice a = d = ±1 is impossible since 
it would lead to more than two independent vectors in the two-dimensional set of fixed points of the 
operators ±5o at t = 0. Thus, we parameterize: a = —d = a, b = j3, and c = 7, where a 2 + /?7 = 1 
and rewrite the relationship ( 15.241 ) in the form: 

02 = 70i - "02 ( 5 - 25 ) 

with the inverse relationship: 

01 = aip! + Pip 2 

02 = 70i - aV>2- ( 5 - 26 ) 

We now wish to solve the linear equation (15.81 . where the right-hand-side vector F(i) is asymptotically 
To-periodic and reversible with respect to <So, i.e. such that 

SoF(-t) = F(i). (5.27) 

By construction, we have 

F = /o0 O + A 00 + 5101 + 5202 

ip = a o o + «i0o + A01 + $202) 
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where the functions gi(t)e A() ' t ' and g2(t)e x °^ are bounded as t — > — oo and asymptotically To- 
periodic. Let M(t) = [4> , 4> ,(f) 1 , (j} 2 ], such that the coordinates of F and ip in the new basis are 

F(t) = M(t)- 1 F(t) = {f ,fi,g 1 ,g 2 ) T (5.28) 
tp{t) = M(t)- 1 <p{t) = (a ,a 1 ,p 1 ,fo) T . (5.29) 

The system (15.81) reduces to the simple system 

—— = -ktai + Jo 
at 

dfr 

-IT = 91 

dp 2 

Since the system (I5.lt is conservative (that is tr (DN(0 o )) = 0), the Liouville's theorem implies that 

det M (t) = det M(0) ^ 0. (5.30) 
Solving the uncoupled system for (3\^{t), we find that 

ft (t) = / 5 2 (s)«!s, 

/9i(t) = f gi (s)ds + (3 w , 
Jo 

such that the functions Pi(t)e~ Xo ^ and /^(i)^ '*' 316 bounded as i — ► — oo and asymptotically To- 
periodic for any constant Pxq. The same property holds for cp(t) as t —>■ — oo. Now, we consider 
the opposite limit t — > +oo. Thanks to the relations (I5.25t . we notice that the functions (agi + 
75 f 2)(i)e _A ° < and (/3<7i — ag2)(t)e Xot are bounded as i — ► oo. Moreover they are asymptotically 
To-periodic if F(t) has this property. Now we can write the function j3\4>i + /?2</>2 as 

/5i0i + /?202 = iPi(t)(J (agi(s) +ig2{s))ds + J jg2(s)ds + af3 10 \ + 

+^2(1) i^J (Pgi(s) - ag 2 (s))ds + J -ag 2 {s)ds + /5/?i ^ 

and it is easy to check the boundedness of the function (3\cj) l + /32</>2 as * provided that 

/ {(3gi(s) - ag 2 (s))ds - ag 2 (s)ds + (3(3 10 = 0. (5.31) 

If /3 7^ 0, this equation determines the value of /3io- If (3 = and a = +1, the value of /?xo i s not 
determined but equation d5.31l > is satisfied identically since 52 (i) is odd, thanks to the identity ( 15.341 ) 
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below. The case = and a = — 1 is excluded by Assumption |2] Thus, the function (3i4> 1 + (3 2 (f) 2 
is bounded as t — ► +00 and asymptotically To-periodic. It remains to show that the reversibility 



constraint (I5.27t implies the anti -reversibility constraint for the solution ip(t): 

S cp(-t) = -cp(t). (5.32) 

By construction, the functions fo(t) and f\(t) are even and odd respectively. Since they are defined by 
a multiplication of an exponentially growing function with an exponentially decaying one of the same 
rate and with periodic factors, these functions are bounded and asymptotically To-periodic. Solving 
the uncoupled system for ao(t) and a\(t), we have 



OLlit) 



«i(0)+ / fi(s)ds, 
Jo 



a o(t) = / [fo(s) - fciai(s)] ds. 





Since fo(t) is even and fi{t) is odd, it is clear that a\{t) is even and ao(t) is odd. Moreover, since 
flit) is odd, it has asymptotically a zero average at infinity. Therefore, the function a\{t) is bounded 
and asymptotically To-periodic at infinity. Now we can find a unique «i(0) such that the even function 
fo(s) — kiai(s) which is asymptotically To-periodic at infinity, has in addition a zero average at 
infinity. Then this gives a bounded function ao(t) which is asymptotically To-periodic at infinity. Now 
for the two last components of F(i) we have 

<&(-t)^i(t) + g 2 (-t)^ 2 (t) - <7i(i)0!(*) " 92(t)<t> 2 (t) = 0, (5.33) 
which corresponds to the identities 

gi(t) - agi(-t) - jg 2 (-t) = 0, 

g2{t) - I3 gi {-t) + ag 2 {-t) = 0. (5.34) 

We need to show that this implies 

Pi{t)+a(3 1 {-t) + 1 f3 2 {-t) = 0, 
fa(t) + m(~t) ~ afo(-t) = 0. 

The identities ( 15.341 ) lead to 

0gi(s) - ag 2 (s) = g 2 {-s) (5.35) 
hence the constant 0iq takes the form 

Pfan = (a-1) ! g 2 (s)ds. (5.36) 

J —00 
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We observe that the identity (I5.36t is satisfied when (3 = and a = 1. Now we check that 

/3i(t) + a/3i(-i)+ 7 /3 2 H) = 



/*oo 

/ {5i ( s ) - agi(-s)}ds + 7 / g 2 (-s)ds + (1 + a)/?i 

7/ g2{~s)ds + (1 + a)/?i = (1 + a)/3i + 7 / g2{s)ds. 
Jo J -oo 



If (3 = and a = 1, the expression above is zero under a special choice of /?io- If (3 / and 
7 = (1 — a 2 )/ (3, it is zero because 



1 + a 



(l + a)/?io + 7/ 92{s)ds = — ^— ( /3/3 10 + (!-«) / ^(s)ds ) = 0. 







Similarly, 



g 2 {s)ds + (3f3 l0 - a g2(s)ds +/ [a^ 2 (-s) - /?gi(-s)] cte 



/3/3io + (l-a) / 52 (s)ds = 0. 
J —00 



Hence the statement is proved in all cases except the case (3 = and a = — 1, which is excluded by 
Assumption |2| ■ 

Using analysis of the truncated vector system (15.11) . we can now prove persistence of bounded solutions 
of the full system. Since there exists a smooth mapping tp(4>,s) in the hyperbolic subspace (see 
Theorem^, we rewrite the system (13.271) in the form, 

= N(0) + ViR(0,V>(<M), (5-37) 

where tf)(4>, e) = O(^). The perturbation term is non-local, which forces us to use a non-geometric 
proof for showing the persistence of heteroclinic connections to periodic solutions near equilibria u± . 
Indeed, let us decompose </> as 

= 0o + V (5-38) 

and define 

Ri(0 o ,v) = N(0 O + (p)- N(0 O ) - DN(cf) )<f = 0(|M| 2 ), (5.39) 

then <p satisfies 

C(t)cp = Ri(0 o , <p) + ^R(0 O + <p, ip(cf> + <f, £)). (5.40) 

Using Theorem 121 and the smooth dependence of R on (cf), tp), we can then use the implicit function 
theorem near <p = and find a unique solution ip 6 C°(M,M 4 ) for small enough e. Moreover, this 
solution (p(t) is a To-periodic function at infinity, which magnitude is bounded by a constant of the 
order of 0{yfe). Using these observations, we assert the following theorem. 
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Theorem 3 Let Assumptions\l\and\2\be satisfied and 4>o{t) b e a solution of the truncated system ( I5.il ). 
which is heteroclinic to small To-periodic solutions u± (t). Then, for small enough e, there exists a 
unique solution <fi(t) of the perturbed system (15.371 1, which is yfe-close to <j>o{t) and is heteroclinic to 
small To-periodic solutions of the system ( 15.371) . 



We note that uniqueness is a remarkable property of the perturbed system (I5.37t . It is supported here 
due to the choice that the function <j>(t) has the same asymptotic period at infinity as the function 

o (t). 

6 Conclusion 

We have derived a scalar normal form equation for bifurcations of heteroclinic orbits in the discrete 
Klein-Gordon equation. Existence of the center manifold and persistence of bounded solutions of the 
normal form equation are proved with rigorous analysis. Bounded solutions may include heteroclinic 
orbits between periodic perturbations at the equilibrium states and true heteroclinic orbits between 
equilibrium states. Our numerical results indicate that no true heteroclinic orbits between equilibrium 
states exist for two important models, the discrete </> 4 and sine-Gordon equations. Double and multiple 
kinks (between several equilibrium states) may exist in the case of the discrete sine-Gordon equation. 

We have studied persistence of heteroclinic anti -reversible connections between periodic solutions near 
the equilibrium states, but we have not addressed persistence of true heteroclinic orbits between the 
equilibrium states. It is intuitively clear that non-existence of a true heteroclinic orbit in the truncated 
normal form must imply the non-existence of such an orbit in the untruncated system at least for 
sufficient small e. Conversely, existence of a true heteroclinic orbit in the truncated system should 
imply a continuation of the family for the full system on the two-parameter plane (c, h). Rigorous 
analysis of such continuations is based on the implicit function argument for a split function, which 
becomes technically complicated since our full system involves the advance-delay operators, where 
the resolvent estimates are not simple. We pose this problem as an open question for further studies. 

A Truncation of the vector normal form from HI95I1 

The eigenvectors for the quadruple zero eigenvalues satisfy the reversibility symmetry relations ( 13.121 . 
The general reversible normal form for this case was derived in II95L We will show here that this 
normal form can be reduced to the fourth-order equation ( 12.91 ) under an appropriate scaling. The 
general reversible normal form is explicitly written in the vector form (see [195 1 for details): 

Ci = & + u 3 P 1 (u 1 ,u 2 ,u< l ), 
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£ 2 = ^,3 + Q3Pl(ui,U 2 ,U 4: ) +U 1 P 2 (u 1 ,U 2 ,U i ) +U 2 P 3 (U2,U 4: ), 

£3 = £4 + r 3 Pi(m,U2,u 4 ) + qiP 2 (ui,u 2 ,u 4: ) + q 2 P 3 {u 2 ,u A ) + n 3 P 4 (ui, n 2 , n 4 ), 

£4 = s 3 P 1 (ui,u 2 ,u 4: ) + riP 2 (ui,U2,u 4: ) + r 2 P 3 (u 2 ,u 4 ) + 93^4(^1, «2, ^4) + u 2, W4), 

where Pi,2,3,4,5 are polynomials in variables: 

ui = £1, u 2 = £ 2 2 _ 2^3, «a = £ 2 3 " 3666 + 3^4, 

w 4 = 3£ 2 2 £ 3 2 - 6£ 2 3 e 4 - 86^1 + i8eie 2 £ 3 e 4 - 9£ 2 £ 4 2 , 

9i = 6, 92 = -3fi& + 6&, 93 = 3£i&& - 4&f| + £& 3j 
ri = 6, r 2 = -3£ 2 £ 4 + 2£ 2 , r 3 = -3£i£ 3 £ 4 + 3£ 2 2 £ 4 - £ 2 £ 3 2 , 

and s 3 = 3£ 2 £ 3 £ 4 — |£| — 3£i£|. If the zero equilibrium point persists in the vector normal form, then 
Ps(0, 0, 0) = 0. The linear part of the vector normal form corresponds to the truncation of Pi 2,3,4,5 as 
follows: 

Pl=0, P 2 =/ii, P 3 = P 4 = 0, P 5 = j U2«l, 

where fii and \i 2 are two parameters of the linear system. The two parameters (/xi, fj, 2 ) must recover 
the Taylor series expansion ( I2.6t of the full dispersion relation ( 12.11 ) in scaled variables. Since the 
dispersion relation from the linear part of the vector normal form is 

A 4 - 3^iA 2 + n\ - n 2 = 0, 

the comparison with the Taylor series expansion ( I2.6t shows that 

w = 4e 7 + 0(e 2 ), fi 2 = e 2 (16 7 2 - 12r) + 0(e 3 ). (A.l) 

The nonlinear part in the dynamical system (I3.14t is scaled as e 2 near the bifurcation point (I3.13t . such 
that 

Pl=£ 2 Pl, P 2 = /Xi+£ 2 P 2 , P 3 ,4 = £ 2 P3,4, P 4 =£ 2 P4, P 5 =^1+£ 2 P 5 , 

where tilded functions Pj depend on the same variables as the original functions Pj. Let £1 = y/e£ 
and use the scaling transformation: 

The scaled version of the vector normal form is written explicitly as follows: 

<f>[ = H + e 3 u 3 Pi, 

4>' 2 = (j) 3 + 47^1 + e<\> x P 2 + e 2 n 2 P3 + e 3 q 3 Pi, 

cf)' 3 = ^ 4 + 4 1 (f )2 + e(l) 2 P 2 + e 2 (q 2 P 3 + u 3 P 4 ) +e 3 r 3 P u 

<Pi = 4 7 3 + (I67 2 - 12r)& + P 5 + e&A + e 2 (r 2 P 3 + g 3 P 4 ) + e 3 s 3 Pi, 
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where tilded variables stand for the scaled versions of the original variables. Let G{4>\) = Pb((/>i , 0, 0). 
The truncated version of the vector normal form (with e = 0) is equivalent to the scalar fourth-order 
equation: 

(f)^ - 12 7 (/>" + 12r^ = G(0), (A.2) 



where = 0i. The comparison with the scalar fourth-order equation i2.9l shows that G(<p) = 
— 12r (F(4>) — cj)), where -F(0) is the nonlinearity of the Klein-Gordon equation d 1 .2b . We note that 
the vector normal form can be simplified near the three curves of bifurcations of co-dimension one 
(see 1195 1 ). The truncated normal form (IA.2t takes into account the full problem of bifurcation of 
co-dimension two. 



B Proof of existence of center manifold in the system (I3.27I) - (I3.28I) 



We shall rewrite the system of coupled equations (I3.27t and (I3.28t in the form: 

</> = N(0) + ViR(0,-0), 
j> = A,o^ + \/iF(0,V), 



(B.l) 
(B.2) 



where <f> G M 4 is a function of (i = yfeC, and xj> G V is a function of The vector fields are given by 



N(0) = (0 2 ,0 3 ,0 4 ,12703 - 12rF(0!)) T , F(0,V>) = gF (p), 



and 



R(0,-0) 



v ^ / 



where g = g — 703 + tF(0i) and g is given below (I3.27t . We shall rewrite g and g in an equivalent 
form, which is used in our analysis. It follows from identities (I3.2t . representation (I3.18t -( 0.19t and 



the scaling (13.26b that 

5±U 3 = 0i ± yF0 2 + \efo ± ^ 3 / 2 4 + e 372 ^, 
5±l7i = 0i(Ci±\^)+£ 3/2 <^i, 

such that the equality = <5±C/i results in the identity, 

03- 2^-01 + ^(<5 + + <T)^ 3 = VF(5+^i+<5_^i-2^i) + 

[0i (Ci + \/e) + <MCi ~ Ve) 
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As a result, the expression for g below (I3.27t can be rewritten as follows: 

g = ^(S+fa + <5-Vi " 2 ^i) + 7e -1 (0i(Ci + \^) + 0i(Ci - - 20i) - r(0 x + e 3/ Vi + Q) 
and 

Q = Q(0i(Ci ~Ve} + e 3/2 5-rp u 0i + £ 3/ Vi, 0i(Ci + VS) + £ 3/2 5+V>i)- 
Using the system (IB.lt . we obtain the following representation: 

= / / (<fe( s ) - —g{s))dsdT, 



such that g is rewritten in an implicit form: 



g = 7 ^(5+^i + (5_V'i-2V'i)-r(</) 1 +e 3 / 2 ^i+Q) 



/ / (03(a) " -g(s))dsdT. (B.3) 



In order to obtain an equivalent expression for g, we use the system (IB.lt and obtain 

-4>3(Ci)+^ 1 / M s )dsdr = e~ l / / / <p 4 (r])dr]dsdT. 

It is also clear that 

Q(MCi-V^),MCi),MCi + ^))-Qo(MCi),MCi),MCi)) 

Cl+VE [dsQ(MCi - >/e),MCi), - d x Q{fa(s - v/i), 0i(Ci), M<i))M* ~ ^)) **■ 

Ci 

As a result, the expression for g = g — 703 + tF(4>i) is rewritten in an equivalent form, 

g = 7Ve (S+ipl + S-ipi - 2-0i) - -7 / / g(s)dsdr 

-r(e 3 / 2 -0i + Q) +ie- 1 f^* [ [ fa^dqdadr, (B.4) 



where 



Q = f Cl+ ^ £ [d 3 Q[M(i-V^),M(i),Ms))Ms) 



Ci 



- d 1 Q (Ms ~ Vi), 0i (Ci), 0i (Ci)) M* - y/e)] ds 
+Q(0i(Ci -y/e) + e 3/2 (5_Vi, 0i + e 3/ Vi, 0i(Ci + + £ 3/2 <5+^; 
5>i (Ci 0i,0i (Ci + V^))- 
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We observe the fundamental property that for 

IMI + lhMI<M 

then there is M > such that 



Ml + 



9 



< M. 



(B.5) 



(B.6) 



In order to prove Theorem ^ we replace the system (IB.2t for VKC) by an integral formula. However 
this is not a simple task, since the estimate on the resolvent of Cx,q is not nice on its third component. 
So, it is necessary here to use the fact that the nonlinear term in the system (IB.lt - dB.2t does not depend 
on ^(Cp), when the system is rewritten in the new formulation (which was precisely the aim of this 
manipulation). For the study of solutions bounded for ( £ 1, we use the two first components of the 
following implicit formula for ip: 

^(C) = -Ve / e^ Cl ' ^-^P + F(ct)(s),tp{s))ds + y/i / e^ --^-^ P^F^s), ip(s))ds, 

Jr J-oo 



'C 

where, by definition 

e (£i,o,+)t p + Y 

e (£i,o,-)* p y 
The curves T± are defined by 



(B.7) 



Lf e A*( A - J C li0 )-ipdA J 
f e At (A — £i o) _1 FdA, 



t < 0, 

t > o, 



r. 



C+ U C + U L+, 

C_ u(Tul_, 



where C± are curves on the complex plane defined by 

C± = I A = ±x + iy : y = ccosh(ax), a > 1, c > \/8, \x\ > £ > o| 



and 



(B.8) 



L± = {A = ±^ + iy : \y\ < ccosh(a£)} . 

The curve F is oriented with increasing Im(A), while the curve r + is oriented with decreasing 
Im(A). It is clear from the proof of Lemma lXTl fhat both curves T± lie in the resolvent set of £i o and 
that on these curves we have (as on the imaginary axis for A = ik, \k\ large enough), 



\Di\ >ci|A| 2 , 



(B.9) 



where D x(A) = 2 (cosh A — 1) — A 2 . We shall complete the system (IB .21) with an explicit formula for 
ip3(C,p). In order to study integrals in the integral equation dB.71 ). we compute the expression 



R(A) = (A-£ li o)- i F (p). 



(B.10) 
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Using solutions of the resolvent equation fl3.16t - d3.17t for F = Fo(p), we obtain the following ex- 
pressions: 

— 34 24 

F((, A) = -- - —(14 + cosh A) + — (cosh A - 1), (B.ll) 



and 



f P l s (l-5s 2 )e A ^ds 
Jo 5 



2p(l - 5p 2 ) 2(15p 2 - 1) 2e Ap Up 



12 



+ - : a h — + — (1 - e Ap ) 

5A 5A 2 ^ 5A 2 A 3 A 4 ; ' 



As a result, the first two components of R( A) satisfy the following bounds for A G T±: 



Ri(A) - 
i? 2 (A) 



5A 2 
3 

" 5A 



< 



< 



Ci 



| A |4-l/a 
C 2 



|A|3-l/a 

for some constants C12 > 0. Since the following integrals cancel for t < (below we consider 
integrals on T + , analogous results hold on T_ for t > 0), 



1 

2iTT 



1 



we obtain 



and 



d 



1 

2ivr 

1 



At 

A 2ivr 



e At Rj(A)dA 



At 

| + ^A = 0, 



< c,- 



-e(a-t) 

a — t 



a 



1 - t 



(B.12) 



(B.13) 



where j = 1,2 and t < min(0,a — 1) (note that a > 1). As a consequence, for a > 1, the two 
first components of the integral ^ J r+ e A *R(A)dA and ^ J r e A *R(A)dA are differentiable real 
functions respectively for t < and for t > 0. We note that no similar estimates can be obtained for 
the third component of R(A) for t < 0. Now using the identity 



A(A - A,o) -1 F = (A - £ li o)- 1 £i )0 F + F 



(B.14) 



one can show that a bounded solution (?/>i, ^2) G Cj 



2 ) of the two first components of the integral 



equation flB.7t is completed by the following explicit formula for ip3(C,p) 



MCp) = ipi{C+P) 
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/ +P (C + P - r) (l - |(C + p - r) 2 ) (£±|±I) dr, (B.15) 



such that -03(C) 0) = V'i(C)- Therefore, for bounded solutions on £ € R, the system (IB .21) with the 
relation dB. 151) reduces to a two-dimensional system. The argument of [ VI92 , p. 145] can then be used. 

Coming back to the system flB.ll) - flB.2t . the implicit linear equation flB.3t for g can be solved for small 
enough e and the system (IB.lt - dB.2l) is hence equivalent to the coupled system d3.27t - d3.28t . From 
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the above study, the two first components of equation iB.lt . completed by (IB.15t have a meaning for 
-0 in C°(R, V), for any known function (ft G C^(]R, IR 4 ). Moreover, for e small enough, and thanks to 
dB.61 ). the implicit function theorem applies, and one can find a unique solution of the integral equation 
( IB .11 for ip £ C°(1R, V), which is bounded by a constant of the order of 0(^/e) under the condition 
that 



where M is an arbitrarily fixed number. Moreover, because of the reversibility of the system (IB.lt - 
dB.2t and oddness property of the function Q(v,u,w), if <fi(Ci) is anti-reversible with respect to So, 
then the unique solution is also anti-reversible with respect to 5. 

It remains to prove the last sentence of Theorem [2 When <p(Ci) is a periodic function, there exists a 
unique periodic solution of the system (IB.2t . since the system (IB.lt - flB.2t is invariant under a 
shift of ( by Tj yfe. When 0(Ci) is an asymptotically periodic function at infinity, the solution VKC) is 
then asymptotic to the corresponding periodic solution of the system dB.21) . This follows from the fact 
that the formula (IB.15t for the third component of can be considered in the limit £ — > ±oo and 
from the two-dimensional equation ( IB.7t for the two first components of tf)(C), for which the solution 
depends continuously on 4>((i), as for an ordinary differential equation. 



The Stokes constants are used in rigorous analysis of non-existence of homoclinic and heteroclinic 
orbits in singularly perturbed systems of differential and difference equations IT00I1 . If the Stokes 
constant is non-zero, the heteroclinic orbit of the unperturbed problem does not persist beyond all 
orders of the perturbation series expansion. We will develop here computations of the Stokes constant 
for the fourth-order equation (I4.lt in the limit a — ► oo, in the context of the 4 model. We will show 
that the Stokes constant is non-zero, which implies that the graph of K(a) on Figure 3(a) remains 
non-zero in the limit a — > oo. 

Let us introduce a small parameter e, such that e = cr _1 , and apply rescaling of t as x = ^Jet. The 
fourth-order equation (I4.lt with F = 0(1 — 4> 2 ) takes the form of the singularly perturbed ODE: 



At e = 0, the ODE (IC.lt corresponds to the stationary problem for the continuous model with the 
heteroclinic orbit in the form <fio{x) = tanh(x/\/2)- There exists a formal solution to the perturbed 
ODE (IC.lt in the form of the regular perturbation series, 




(B.16) 



C Stokes constant for the fourth-order ODE (14.11) 




(C.l) 




k>l 



(C2) 
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where an odd function 4>k (x) is uniquely defined from the solution of the linear inhomogeneous prob- 
lem: 

{-d 2 x + 2 - 3sech 2 -^Q <f) k = i?fc(0 fc -i,0fc-2, -Ao), k>l, (C.3) 

where H k are correction terms which are odd exponentially decaying functions of x £ R. Each term 
of 4>k{x) is real analytic function of x G R and is extended meromorphically to the complex plane, 
with pole singularities at the points x = ?Qk™te. t n £ Z. The Stokes constants are introduced after 
the rescaling of the perturbed ODE dC.lt near the pole singularities and replacing the formal solution 
dC.2l) with the inverse power series. We rescale the solution as follows: 

iri .. . if>(z) 

X=—= + €Z, 4>(X) = , 

y/2 e 
where ip(z) satisfies the regularly perturbed ODE: 

^(iv) + yll _ ^3 + = Q (C _ 4) 

We are looking at the formal solution of the ODE dC.4t with the perturbation series 

i,{z)=Mz) + ^ k Mz), (C.5) 
fe>i 

where each term ifj k {z) is the inverse power series of z, starting with the zero-order solution: 

Matching conditions between dC.2t . dC.5t . and dC.6t imply that a\ = v2 and 02 = 0. The main 
result of the beyond-all-order perturbation theory is that the formal series dC.2l ) are dC.5l > diverges if 
the inverse power series dC.61) diverges. Divergence of the formal series implies in its turn that the 
heteroclinic orbit ^o(^) = tanh(x/\/2) is destroyed by the singular perturbation of the ODE dC.lt . 
Divergence of dC.6t is defined by the asymptotical behavior of the coefficients a m as m —> 00 (see 
ITP051 for details). 

In order to show that the inverse power series (IC.6t diverges, we find the recurrence relation between 
coefficients in the set {a m }^ =1 : 

m(m + l)(m + 2)(m + 3)a m + (m + 2)(m + 3)a m+2 - aia k a m+A _i_ k = 0, (C.7) 

Z+fc<m+3 

where m > 1, ai = v2 and 02 is arbitrary. Due to the translational invariance, we can always fix 
02 = 0. Setting a m = (— l) n (2n)\b n for m = 2n + 1, n > and a m = for m = 2n, n > 1, we 
reduce the recurrence equation dC.7t to the diagonal form: 



b n+ i = b n + Yl (2n + 4)! L kb k b n+1 ^ k , n > 0, (C.8) 

l+k<n+l ^ 71 '' 
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where b = \/2. Since 

6 

1 - - r > 0, n > 0, 

(2n + 4)(2n + 3) ~ ' 

we prove that the sequence {b n }^ =0 is increasing, such that b n+ \ > b n > bo = y/2. Therefore, the 
sequence {6 n }^L is bounded from below by a positive constant, so that the Stokes constant in the 
beyond-all-order theory is strictly positive. 

Similar computations can be developed for the fourth-order equation (14.11) with F = sin(^). They are 
left for a reader's exercise. 



D The normal form for the inverse method of [FZK99] 

It was shown in [FZK99| that the nonlinearity of the discrete Klein-Gordon equation (11.11) can be 
chosen in such way that it yields an exact travelling kink solution u n (t) = 4>{z), z = n — ct for a 
particular value of velocity c = s. In application to the discrete (ft 4 lattice, one can look for the exact 
solution in the form (f>(z) = tanhQuar), where /j, is arbitrary parameter, and derive the explicit form 
of the nonlinearity function / = f(u n ) parameterized by s and [i. (See eq. (27) and refs. [13,14] in 
[FZK99].) Since we are using a particular form of the nonlinearity f(u n -i,u n ,u n+ i) in the starting 
equation (ll.lt . we transform the nonlinearity from [FZK99| to the form: 

f( Un ) = Un (l-u 2 n )\±^, (D.l) 

where 

s 2 fi 2 tanh 2 [i 2 

a= l — Z2 rT> P = tanh n, 

tanh \x — s L n z 

and h 2 is related to the parameters \i and s by means of the relation: 

h 2 = 2 (tanh 2 jj, - s V) . (D.2) 

We note that the nonlinearity (ID. II ) violates the assumption (ll.6t of our paper, since it depends implic- 
itly on the lattice parameter h 2 . 

When s = (a = 0, f3 = \h 2 ), the nonlinearity (ID.lt represents one of the four exceptional non- 
linearities with translationally invariant stationary kinks (see [BOP05|). We consider here the limit 
s — > 1 from the point of normal form analysis. We will show that existence of the exact travelling kink 



solutions 4>(z) = tanh(^z) in the differential advance-delay equation (ll.lOt is preserved within the 
reduction to the scalar fourth-order equation \2.% . Let 

s 2 = 1 + e~f s , c 2 = 1 + 67, h 2 = e 2 r. 



The parameter fj, is defined from the transcendental equation dD.2t in the form, 

H = ^fen s + 0(e), 
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where /j, s is a real positive root of bi-quadratic equation: 

-l4 + 2 ls £ + r = 0. 

Two real positive roots fi s exist in the domain < r < §7^ for 7^ < and no real positive roots exist 
beyond this domain for r > 0. The nonlinearity (ID.lt is reduced in the limit of small e to the form, 

F{u) = u(l - u 2 )(l + ^u 2 ) + 0(e), 

r 

such that the normal form ( I2.9t becomes now 

l^v) _ # // + T0(1 _ ^ + 2m 4^3 (i _ 02) = Q (D 3) 

It is easy to verify that the scalar fourth-order equation dD.3t has the exact heteroclinic orbit 4>{z) = 
tanh(/i s z) for 7 = 7 S and r > 0. Since there are two roots for /j, s , two heteroclinic orbits <f)(z) exist. 
These two orbits disappear at the saddle-node bifurcation as r — > jj 2 . Within the numerical method 
of our paper, the heteroclinic orbits can be detected as isolated zeros of the split function K(a). We 
conclude that modifications of assumptions on the nonlinearity of the discrete Klein-Gordon equation 
( ll.lt result in modifications of the scalar fourth-order normal form d2.9t . which may hence admit 
isolated heteroclinic orbits. 
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